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entire Brillouin zone (and not with a low-energy effective action) we reduce 



I successively the cutoff A about the Fermi surface and follow the renormaliza- 

(N 

, tion of the coupling U as a function of three energy-momenta. We calculate 

I the intrinsic scale Tco where the renormalization group flow crosses over from 

^ I the regime (A > Tco) where the electron-electron (e-e) and electron-hole (e- 

' h) terms are equally important to the regime (A < Tco) where only the e-e 

Q I term plays a role. In the low energy regime only the pairing interaction V is 



marginally relevant, containing contributions from all renormalization group 



I steps of the regime A > Tco- After diagonalization of Va=Tco^ identify its 

most attractive eigenvalue Amin. At low filling, Amin corresponds to the B2 
representation (d^y symmetry), while near half filling the strongest attraction 
occurs in the Bi representation {d^2_y2 symmetry). In the direction of the van 
Hove singularities, the order parameter shows peaks with increasing strength 
as one approaches half filling. Using the form of pairing and the structure 
of the renormalization group equations in the low energy regime, we give our 
interpretation of ARPES experiments trying to determine the symmetry of 
the order parameter in the Bi2212 high-Tc compound. 
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I. INTRODUCTION 



An immense number of recent experiments on high superconductors aims at a de- 
termination of the form of the BCS gap function in momentum space.§ Josephson junction 
experiment sjm measurements of the London penetration depth,S@ and of the Cu NMR re- 
laxation ratei are consistent with a (i^2_j^2-gap. Particularly interesting are the ARPES 
data which provide rather precise information about the detailed angular dependence of the 
amplitude of the gap function. The experiments on the Bi2212 compound show that the 
order parameter is maximal along the (0, vr) directionIBi and that its amplitude in the (vr, tt) 
direction seems to attain a nonzero value at a new critical temperature below Tc.i 

The d symmetry of the gap function is generally considered to be a sign of a pairing 
interaction of electronic origin, implying the absence of the standard phononic mechanism 
for superconductivity. The idea of a superconducting state induced by fluctuations of purely 
electronic origin in systems of electrons with Coulomb repulsion is originally due to Kohn 
and Luttingei0 for the case of the three dimensional electron gas. Similar effects exists 
in a two-dimensional electron gas, and generally they depend strongly on the form of the 
Fermi surface. Perturbative calculations of the four point vertex for a weakly filled band 
in the Hubbard model show that the model is instable against d^y superconductivityS. 
Quantum Monte Carlo calculations on the same model, in the vicinity of metal-insulator 
transition, show that the attractive pairing interaction of d^i-y^ symmetry is dominant ,0, in 
agreement with the earlier arguments that antiferromagnetic fluctuations are the mediator 
of pairing interactionsEHli. Direct evidence for dx2_y2 superconducting ordering has however 
not yet been found in quantum Monte Carlo studies. Antiferromagnetic fluctuations become 
stronger and stronger as one approaches half-filling. For weakly interacting electrons, these 
fluctuations are associated with the In^ T-divergence of the electron-hole (e-h) loop diagram, 
caused by the nesting property of the Fermi surface and the van Hove singularities. On the 
other hand the BCS fluctuations, characterized by the electron-electron (e-e) loop, normally 
only linear in logarithms, crosses over to a In^ T form in the vicinity of half-filling. The 
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perturbative treatment of an interacting system of electrons should thus be based on the 
summation of all iterations of these two types of loops. The renormalization group is one way 
to do this. Applied to interactions between electrons placed at the van Hove points it gives 
an antiferromagnetic instability at half-filling and superconductivity of dx2_y2 symmetry if 
the deviation of the chemical potential /i from its value at half-filling becomes of the order of 
critical temperature of the antiferromagnetic state.0 The equivalent parquet approach has 
been used for half-filling and also finds the antiferromagnetic state.S'il A direct calculation of 
the zero temperature free energyil up to the second order in the bare interaction Uq confirms 
that at half-filling the antiferromagnetic order is stable, but finds no finite superconducting 
order parameter at any filling. 

In the present analysis we search to know whether the Hubbard model with repulsive 
on-site interaction can lead to superconductivity, and if it does, to what form of the gap 
function. In particular, we are interested in the dependence of the results on the density of 
electrons. It is hoped that the results can help us to clarify the origins of the existence of a 
highly anisotropic BCS gap function in the cuprates. 

The renormalization group technique for fermionic systems in two and three dimen- 
sions has recently been developed, but with very drastic limitations. The Wilsonian mode- 
elimination technique was applied by Shankaiii only to systems with an either isotropic or 
open, perfectly nested Fermi surface. Weinberg!! has written the flow equations for a gen- 
eral case of an anisotropic Fermi surface, but taking into account only the electron-electron 
channel of the flow. Moreover, a common tendency is to do the renormalization group proce- 
dure only for a thin ring of degrees of freedom around the Fermi surface and to linearize the 
spectrum in the radial direction, taking as a starting model the low-energy effective action. 
This makes it difficult or impossible to make statements about the phase diagram of lattice 
models like those relevant for the description of the cuprates, due to the absence of a proper 
description of high-energy degrees of freedom, the elimination of which may considerably 
affect the effective low-energy action. 

We formulate the renormalization group starting from the Hubbard model, with the 



whole Brillouin zone involved in renormalization. We do no tree level scaling before the 
energy shell ±A of available states around the Fermi surface contains no van Hove singularity. 
Once A has become the smallest energy scale, performing a tree level analysis and deriving 
the renormalization group equations, we show that the effective coupling function contains 
some relevant contributions with the origins in the regime with higher A, i.e. from the 
electronic degrees of freedom not included in the low energy effective action. We must add 
that our renormalization group, since perturbative in interaction, can provide uniquely the 
information whether the Fermi liquid is a fixed point or not. If the coupling flows to strong 
coupling, we can say in which direction it flows, for example in the rf-type superconducting 
direction, but we can not say whether another fixed point with finite superconducting order 
parameter with d symmetry exists or not. This kind of problem is well known e.g. from 
the renormalization group in quasi-one-dimensional compounds, where the most divergent 
flow in some direction is always associated with the corresponding long-range order (LRO) 
because already infinitesimal interchain coupling suffices for its stabilization.ii Similarly, the 
dimensionality reason for the non-existence of LRO in two dimensions at finite temperature 
can be ignored as soon as small hopping in the third direction exists. This, however, doesn't 
mean that LRO and Fermi liquid fixed points are the only possible, on this question the 
one-loop renormalization group simply can not give an answer. 

The problem is formulated in section 2 in terms of the effective action and the corre- 
sponding renormalization group flow equation for the coupling function. In section 3, after 
calculation of the crossover energy to the purely electron-electron (e-e) part of the flow, 
we derive the renormalization group equation for the pairing function V. In section 4. we 
diagonalize the pairing V for the case of a non-divergent electron-hole (e-h) channel and de- 
termine the most attractive eigenfunction and the resulting critical temperature as functions 
of the chemical potential /i. In section 5 we give our picture of the gap viewed by ARPES 
experiments. The conclusions are given in section 6. 
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II. MODEL AND FORMULATION OF THE RENORMALIZATION GROUP 

THEORY 



The Hubbard model for a two-dimensional system of electrons on a square lattice is given 

by 

^ = E^k«lk«<xk+ ^f/oE E «L,kl+k2-k3«-^k2 4k3«^kl , (2.1) 

crk o" ki,k2,k3 

where = — 2i(:(cos + cos ky) — fi and the momenta are within the first Brillouin zone. 
For this model we introduce now the action in terms of fermion coherent states represented 
by Grassmann variableJ^ ^aK and ^^k, where K = (k,uj). We will write it in the form 

Jo 2tt 

+ ^E/fn^) E wLtlk3,k,^.i.3^.'i..^.'i..^.i.. , (2.2) 

era' \i=l / ki,k2,k3 

where 0kl!k2,k3,k4 = 11^=1 ©(A — |^kil) constrains all four momenta to run within the energy 
shell ±Ao = 8t around the Fermi surface. The energy and momentum are conserved so that 
K4{Ki, K2, K3) = (tui + LJ2 — ki + k2 — ks). Note that the size of the cutoff is equal 
to the bandwidth, i.e. the whole Brillouin zone is available for integration. Thus, the 
functions have no meaning yet: they become important when the cutoff, reduced by the 
renormalization group, become lower than a tilt from the Fermi level to the band boundary. 
Note that for a non-half filled band the effective phase space is not particle-hole symmetric. 

The renormalization group transformation that we will use, known as the field theory 
approach, is defined as the mapping 

S{Ao,Cl Uo} ^S' = S{Ao ^ Aoe-',ek ^ a, Uo ^ UiK,, K2, K,)}, (2.3) 

where and U depend on / in a way so that the physical properties of S' and 5* are the same 
for energies lower than A = Aoe'K This requirement is fulfilled if all one-particle irreducible 
vertices are invariant under reduction of the cutoff from Aq to Aqc^K The renormalization 
group can be thought of as a set of successive, infinitesimally small steps dl. This allows 
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us to formalize the renormalization group requirement in a set of equations diFi — 0, where 
i = 2,4,6, .... Up to second order in [/, it suffices to consider only r2 and because 
all higher vertices are of higher order in U. The solutions of the equations diT2 = and 
dir2 — give us the renormahzation group flow for and U{Ki, K2, K^). 

Conservation of spin allows us to write the interaction part of the action as a sum of the 
singlet (|a + cr'l = 0) and triplet + cr'| = -s/2) parts: 

s(K4, K^)U\K^, K2, Ks)s{K2, Ki) + i^{K^, K^)U^{K^, K2, Ks)t^{K2, K,), (2.4) 

where s and are the variables of annihilation of the singlet and triplet states 

S{K2,K^) = -L^a*,^,*_,;^,, (2.5) 

to{K2,K^) = -L^^^^^^_^^^ . t±^{K2,K,) = (2.6) 

The singlet state is symmetric and the triplet antisymmetric under exchange of the momenta 
of two particles. Correspondingly, the couphng function U^{Ki, K2, Ks) can be taken to 
be symmetric and U^{Ki, K2, K^) to be antisymmetric under the momentum exchange 
operation X, defined as 

XJ^{K,, K2, K^) = T{K2, i^i, K^) , (2.7) 

T being a function of four energy-momenta which conserves energy and momentum. If T 
possesses time-reversal symmetry 

^(i^i, K2, Ks) = TJ^{K,, K2, Ks) = ^(i^3, K^{Ki, K2, Ks), K,), (2.8) 

which certainly is a property of the vertex, then it is equivalent whether X exchanges Ki 
and K2 or K^i and K^, i.e. T{K2,K^,K^) = J^{Ki, K2, K^iK^, K2, K3)). Formally, C/^ and 
are given by 

U^^^{1-X)U, U'^^{1 + X)U. (2.9) 

On the other hand, the interaction can also be written as a sum of one term with equal 
(cr = a') and one with opposite {a — —cr') spin quantum numbers, with corresponding 



coupling functions named [/y {Ki, K2, K3) and Uj_{Ki, K2, K3), respectively. From two equal- 
spin electrons one can build only a triplet state, which make us conclude that 

U\\ = f/^, (2.10) 

while 

U^ = U = U^ + , (2.11) 

containing the singlet and the triplet interactions. 

Another way to write the interaction is in terms of a charge and a spin part 

U,{Ku K2, K^)C{K2, K,)C{Ks, K^) + U„{K^, K2, Ks)S{K2, K,) ■ S{Ks, K^) , (2.12) 

where C and Si are 

C{K,, K,) = J2 ^.K.'^.K, ; S,{K,, K,) = J2 ^aK-Ao'-^^'K, ■ (2.13) 

The charge and spin coupling functions are 

U, = h^2-X)U, U^ = -^U. (2.14) 

In first order of perturbation theory, these functions determine the renormalization of the 
charge-charge and the spin-spin correlation functions. 

We proceed now with the derivation of renormalization group equations. For simplic- 
ity, we will ignore the renormalization flow of which follows from conservation of 
renormalizing the form of the Fermi surface, the effective mass, etc. This approximation is 
justified in the case of the circular Fermi surfaced. In the anisotropic case, the diagrams 
for r2 have a dependence on the direction of k. Moreover, even a small renormalization of 
the Fermi energy can give important changes of the form of the Fermi surface if one is close 
to half-filling, because of van Hove singularities. For filling not too close to one-half we can 
expect that the essential of the physics is given by just the renormalization of the coupling 
U using the bare dispersion relation which we will call .^k from now on. 
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The Feynman diagrams for r2 and T4± = T4 are given in Fig. 1. The first loop in the 
expression for r4 is of the electron-electron (e-e) and all others of the electron-hole (e-h) 
type. Making use of the relations ( |2.10| ), ( |2.11| ), and (|2^), we get the expression for r4 in 



terms of U and XU. If we write the integration measure of the loop diagrams in the form 

^ ^e/^^, 2.15 
2n J-A J f (s, 4) 

s being the curves of constant energy ^, then dT^ corresponds to the integration of the two 
energy shells of width \A\dl at ,^ = ±A. We obtain the following flow equation 
dU 

— = (3ee{U, U} + MU, U} , (2.16) 

with 

MU,U} = 2f3eh{U,U} - (3eh{U,XU} - f3eh{XU,U} - X[3eh{XU,XU} (2.17) 

The functionals l3ee{Ui, U2} and l3eh{Ui, U2} are the partial derivatives with respect to I of 
the e-e and e-h loops and both are bilinear forms in Ui{Ki,K2,K3) and U2{Ki, K2, K3). 
They read 

(3ee{Ui, U2} = (H{f/i, U2} + E{XUu XU2}) (2.18) 

and 

MUi, U2} = (n{?7i, U2} + ru{Ui, U2}) (2.19) 

with 

-A ^ r ds„ ^ , . _ r+°° duo 1 1 



{2Tiy ^^_J Vy " 7-00 2 7r itU - Z/A Z(-Lc; + tUee) - ^k.-qee 

X f/i(iri, K2, K(^u))U2{K3, K4, i^H) , (2.20) 

1 



n{^i, U2} = j^,T. / (A - |e,.^,,j) r ^ -1 



■X 



Z/A i{uJ + UJeh) -^k.+qe„ 

X f/i(iri, Ks)U2{K4, Kh, K2) . (2.21) 
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The index i/ = +, — symbolizes two energy shells at +A and —A; Vi, stands for ^ = i^A); 
Ue.e = o;! + U2] LUgh = coi - CO3; Qee = ki + ks; Qefe = ki - kg; = {K,uj), where K is 
the momentum running along the path s^. Kee and Kgh are non-analytic functions of the 
momenta, given by 



for Qee = 

1 otherwise 



for qeh = (±7r, ±7r) 

1 otherwise 



Their origin is in the derivatives over A of the products like 

e(A-|ek+cil)e(A-|ek|) 

when = Ck+q- 

III. THE TEMPERATURE SCALE Tco 

A particularity of the renormalization group approach treating e-e and e-h fluctuations in 
more than one dimension is the absence of self-similarity of the problem. In fact, there is an 
intrinsic energy scale which is a function of the band filling. It is associated with charge and 
spin fluctuations coming from the e-h term. We will proceed by estimating the characteristic 
energy scales which appear in (3ee and (3eh, when all four particles are at the Fermi surface, 
with zero energy. If we are exactly at half-flUing, it is known that in the limit a; — > both 
e-e and e-h loops scale like P, which corresponds to the square- logarithmic divergence in 
both channels. This gives an exphcit ~ / dependence in the /?-functionals. Let us suppose 
now that the filling is slightly lower than one-half, i.e. that /i is small and negative. We 
expect two regimes. One is for I ^ ~ ln|8i///|, where the flow is still unaffected by 
the small changes of the Fermi surface due to nonzero n and remains proportional to /. In 
the second regime, where I ^ Ix, the e-e flow is just a constant (i.e. only a In-divergence) , 
while the e-h flow decays exponentially due to disappearance of nesting. Even far from 
half-filling it is possible to define a crossover l^, beyond which the flow in the e-h channel 
disappears exponentially. We can summarize saying that for any filling, is a crossover from 
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a regime where both e-e and e-h loops contribute to a regime where (3eh starts to behave 
hke I3eh ~ A^'^'^^ Here //(q) is positive for all values of the momentum transfer q = qe^. 

To estimate the dependence of Ix on the filling (ra), we consider the static limit of the 
partially integrated e-h loop 



with the momentum transfer q equal to 2k^ in the direction (7r,7r). Note that the energy 
integration is performed over 8t < < 8texp(— /). The derivative of Peh(^q, = 0) with 
respect to / gives the explicit /-dependence in the [3ee functional. Fig. 2 shows Peh{l) and 
diPeh{l) for two different values of /i. It is reasonable to define Ix as the point where diPf.h{l) 
starts to decrease. In the exponential regime the function diPehi^) decays like exp(— //2) 
(i.e. ?7(2k^) = 1/2, valid for any orientation of k^), while the regime I < remembers the 
In^ divergence of Peh at half filling. If we consider PehiO for some large momentum transfer 
different from 2ki7 (giving intersection rather than touching of the initial- and final-state 
Fermi surface) we get a shape like exp(— /) (i.e. ri{q ^ 2ki7) = 1). The dependence lx{{n)) 
is shown in Fig. 3. Near (n) = 1 there is a divergence of the form lx{{n)) ~ In |8t//i((n))| 
because of nesting, while the increasing as the filling goes to zero mirrors the fact that, 
for low density, the Fermi energy appears as the new scale instead of the band-width being 
used. The insert shows the function (n)(/i). 

Once in the exponential regime, f3eh can be neglected after it becomes smaller than 
of its value a.t I = Ix- Putting r] = 1/2, this defines the crossover 



corresponding to the crossover temperature Tco = 8texp(— /co)- Suppose that we now inte- 



l3ee{U} remains in the flow equations, and one has a partially renormalized U{1 = Ico) as 
initial condition. Now we use the fact that Tco/8t is a small parameter, i.e. the inequality 
I'Ckl < Tco determines a thin ring of degrees of freedom, containing no van Hove points, as 




(3.1) 



^co(Ai) = IxifJ') + 2 , 



(3.2) 



grate the flow equation ( 2.16|) from / = to / = Ico- Once / has reached Ico, only the term 
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one can conclude looking at Fig.3. This allows us to rescale the momenta k± = n(k — kp), 
where n is the unit vector normal to Fermi surface, dependent on direction of k. To clarify 
the reason for which a tree-level scaling is not allowed for energies higher than the deviation 
of the Fermi level from van Hove singularity, let us write the phase space integration measure 
in terms of energy (^) and polar angle (6) variables 



with J{^,6) = k{^,6)/v{^,6), k being the radial wave number and v the group-velocity. 
Zeroth order (tree level) scaling tells us via power counting argument to consider J{C,, 6) 
as function of 9, neglecting any ^ dependence about ^ = 0, which is possible if J(^, 9) is 
an analytic function of ^ at the whole shell ±A which, consequently, should contain no 
singularity. 

We can also rescale the frequencies if U{lco) is an analytic function of u in the interval 
±Tco about the Fermi surface, which we assume to be the caseS. In the scope of this tree- 
level scaling, as it has already been shown by Shankarii, the slope of the electronic dispersion 
around the Fermi surface is irrelevant and the two marginal interactions correspond to two 
different constraints on the four-momenta in U. Since any k±- and cj- dependence in U 
is irrelevant, both marginal interactions depend only on coordinates of the zero frequency 
particles placed at the Fermi surface. For the first, "Fermi liquid" or forward interaction, 
the momenta satisfy the equation ki = ks, where the meaning of momenta can be seen from 
the equation (p.2|) . This interaction is slightly (~ Uq) renormalized by the high-energy- 
modes (/ < Ico), and is not involved in further renormalization. The second interaction is 
the pairing potential V, where the momenta satisfy the condition ki = — k2. The pairing V 
depends only on angular coordinates of annihilated and created pairs. Keeping in mind the 
above remarks, we can write the action for the electrons in the ring ±Tco around the Fermi 
surface as 




(3.3) 




ds 



^,{e,9){dr-e)^^{e,9) + 



2nv{9) 
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+ ^S/^/^/^A._«,')V',.,JM')A.,_W + FL}, (3.4) 

where e = k^v{6), closed loop integrations are over the the Fermi surface and FL stands for 
the effective Fermi liquid interaction. Aq-/ o- q^^(^^) is the energy-integrated number of pairs 
defined as 

A^^^^^^^je)^ I ^vl>^^(k)vl/.,(-k + q,e)e(T,„-|(-k + qe,)-nt;(^)|). (3.5) 

Note that the integration measure over small momentum qee goes to zero as / — oo. A form 
similar to ( p.4|) has been used by Weinbergil, but taking Vi=i^^ {9, 6') phenomenologically 
and not as the partially renormalized pairing interaction which we get from Ui=i^^{6i, 02, Os) 
putting incoming particles 1 and 2 to ^ and 6 + n and outgoing 3 and 4 to 6' and 6' + vr. 
Note that the loop integration over s can be understood as the scalar product over "vector 
components" of a "spin", where the number of components corresponds to St/Tcdl'il. 
The integrations have the weight factor l/v{6), v{6) being the anisotropic Fermi velocity, 
what suggests to introduce a new angular coordinate 



rs ds 
J v(e) 

27rNF 



m = (3.6) 



where Np is the density of states at the Fermi level. The function z{6) is shown in Fig.4 for 
few different values of fi. Starting from the new Hamiltonian ( p.4|) we can calculate now the 
function Pee{V} in ^-space and obtain the flow-equation 

diViz, z') = -^i dz"V{z, z")V{z", z') , (3.7) 

where the coordinate z appears instead of 9{z). For initial condition we take 

Vi=i^XO{z),e{z')). 

IV. DIAGONALIZATION OF THE PAIRING POTENTIAL 



To make the differential equation ( p.7|) solvable one has to diagonalize the pairing po- 
tential V{z, 2;'). Hi! Since it is invariant under all symmetry elements of the point group, 
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its most general form in z-space can be written as 

V{z,z')=J2T.V2,nflni^,^'), (4.1) 

7 m,n 

where = {m'y\V\n'y) and fm,n{z^^z') = {m'y\z){n'y\z'). The function {mj\z) is the 

m^^ basis state of the 7-representation of the point group D4^. It is proportional to the 
function cos4mz, sin4mz, cos(4m + 2)z, sin(4m + 2)z, and [cos(2m + l)z ± sin(2m + 1)2;], 



for J = Ai, A2, Bi, B2, and E respectively. Using the flow equation (|3.7|) becomes 



Nj? 

dV"' = -y^V^ (4.2) 

with the initial condition 

yiAl = Ico) = I dzdz'ajz, z')ViJz, z'). (4.3) 



To solve exactly the equation ( |4.2|) one has to diagonalize five infinite dimensional matrices 
V^„(/ = Ico), thus decoupling completely the flow ( |4.2| ) into a set of differential equations 
whose solution is 

Here A labels the eigenvalues within the representation 7. If V^{lco) is negative, the denom- 
inator has a zero at / = /c(7, A) and an instability occurs. 

The renormalization group calculation of the pairing interaction V{lco) is an extremely 
difficult problem because of the interplay between the cutoff A and the geometry of the 
two-dimensional shells dA. However for some range of (n) and Uq (essentially small Uq and 
(n) not too close to half-filling) the problem can be reduced to only one simple integration, 
i.e. 

V{ho) ~ f/o + /o (4.5) 

with 

Io = U^{Pee{lco)+Peh{lco)), (4.6) 
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where Pee(0 is the partially integrated e-e loop, defined in a same way as Pehi}) in eq. (|3.1| ). 
Peh{^) comes just from the last diagram in the expression for r4 in Fig. 1, because the other 
three e-h loops cancel exactly for any U symmetric under exchange. The approximation 



( [4.5| ) is allowed if Iq/Uq is a small parameter. A good estimate of the magnitude of this 
parameter is given by the Ai-part of UoPeh{ico), giving the criterion for the validity of our 
approach in ((n), t/o)-space. The line Iq{^)/Uq ~ 1 is shown in Fig. 5. Below the line the 
approximation is justified. 

By definition of /co, Pehilco) ~ Pehil = oc). Moreover, Peeilco) has no dependence on z 
and z\ since it depends on external momenta only through ki + k2, which we put to zero. 
Consequently, its only nonzero component is (0y4i|Pee|0^i). Thus, for the calculation of all 
other components of VJ^^ we use just the bubble PehQ- oo) with the momentum transfer 
q = \s~f{z) — \i~F{z' + vr). Fig.6(a) shows Peh{l oo) as a function of z and z' for chemical 
potential /i/4t = —0.2. 

The minimal eigenvalues of V{lco) in all five channels, named A^jj^, are shown in Fig. 6(b) 
as functions of /i. These curves indicate which kind of superconducting symmetry becomes 
critical at some given /i. The eigenvalues for each channel are calculated taking only the first 
four harmonics for Ai, A2, Bi and -B2, and the first six harmonics for the E representation. 
The corresponding eigenvectors determine the Fourier spectrum of the gap function. A very 
important result is that the relevant harmonic of the superconducting fiuctuations in the 
Bi channel occurs very close to just cos(22;), being thus determined only by the structure of 
the Fermi surface and not by the interaction. Fig. 7(a) show the instable order parameters 
^^ = —0.5, for /i = —0.31 and A^^ for /i = —0.001 as a function of the Fermi 
surface angle 9. The evolution of the function cos22;(^) (i.e. the first harmonic of -Bi) with 
log(— /i), given only by the dependence of z on 9, is shown in Fig. 7(b). The strength of the 
peaks near the van Hove points increases and the magnitude in the area between the peaks 
decreases with log(— /i). 
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The critical temperature is given by a cutoff for which the most attractive diagonal 
component of V diverges, i.e. 

T, = 8texp[-/e(7,Amin)], (4.7) 

where Amin = min{A^jn}. Fig. 8. shows Ic as a function of log(;u). The critical temperature 
decreases extremely fast as we go away from the half-filling. An increase of Uq could save 
the situation, but in that case our perturbative method ceases to be sufficient (see Fig. 
5). Since the cuprates are superconductors for fillings quite far from one electron per site 
((n) ~ 1 — 0.17), this result means that the small-t/ Hubbard model cannot describe these 
systems quantitatively. However, the model gives very precious informations about the form 
of the gap function in the i?i-instable regime, which will not change considerably with 
increasing Uq, as long as cos(2z(6')) is the dominant attractive harmonic in V{lco)- 



V. MIXED SYMMETRY SUPERCONDUCTIVITY 

Once the renormalization flow has been integrated for I < Ico, assuming that the interac- 
tion did not diverge earlier in the antiferromagnetic channel, the detailed angular dependence 
of the superconducting gap function can be easily found. In general, a superconducting state 
with the symmetry corresponding to the lowest of the eigenvalues A^j^ will be formed. How- 
ever, when two of the A^j^ are close to each other, a more complicated situation can occur: 
for definiteness, consider the region 0.206 < \fi/4:t\ < 0.276 in Fig. 6(b), where the B2 eigen- 
value is the most attractive after the Bi. Let us suppose that Bi order of the simplest form 
A^^ ~ cos2z has formed and that the temperature is close to T^. Among the remaining 
symmetry channels, B2 is the only one which can give a large gap function in the node-points 
of A^^ and zero in the points where A^^ is maximal. Consequently, we expect that the flow 
of the type ( [4.4| ) with 7 = i?2 will not be strongly affected by the existing Bi order. Consid- 



ering to a first approximation the two flow equations (7 = Bi and 7 = B2) as independent, 
and taking only the first harmonics of the Bi and B2 representations, we can construct the 
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relevant part of the pairing interaction which gives two phase transitions, one with Bi and 
the other with i?2-symmetry: 

Vi^^{9i,92) = V'^^''>- cos2z{9i) cos2z{92) + V^^^'>- sm2z{9i) sm2z{e2), (5.1) 

TT TT 

where all details of the Fermi surface are contained in the dependence of z on 9. From 
eq.( [4.4| ) one finds that the ratio between two critical temperatures is given by T'^/T^. = 
exp[-27r(l/1/(-^i) - 1/V^^''^)/Nf]. From Fig.6(b) note that the ratio T^/T^ is very sensitive 
to the variation of the chemical potential. The gap function resulting from ( |5.1D has the 
form 

A{9) = As.e^'^' cos 2z{9) + Ab^c''^' sin 2z{9) , (5.2) 

where A^^, Ab2, 0i and 02 are real. These parameters can be determined minimizing the 
mean-field expression for the free energy per siteS 

F = -^E0(^co - la|)lncosh^ + |Ai|Vv(^^) + \A2\VV^^^\ (5.3) 



where E'k = Y^k+ ^^"^ theta function constrains the momentum summation 

to run only over the states within the energy shell ±Tco about the Fermi surface. The 
minimization of F with respect to cos(0i — 02) gives 

= 01-02 = ±^, (5.4) 

i.e. the resulting gap function is of the type Bi ± ii?2. It is interesting to remark that the 
same kind of gap function has been obtained by Laughlin using the anyon picture.i3 The 
particularity of this gap function (and of any gap consisting of two different symmetry terms 
with a phase difference of ±vr/2) is that it breaks time reversal symmetry. 

We can now try to understand recent ARPES measurements by Jian Ma and coworker^ 
on the Bi2212 compound. From their experiment it appears that two superconducting 
instabilities occur; the first one is at T = Tc and has probably the Bi symmetry. The second 
instability occurs at = O.SlTc, it introduces a nonzero gap at the points 9 = (2n+ l)vr/4, 
i.e. halfway between the corners of the half-filled Fermi. The function measured by ARPES 

17 



is |A(^)| and in the picture discussed above has no zeros and minima on the diagonals of 
the Brillouin zone if both A^^ and are finite, |AbJ > jA^jl and (f) = ±7r/2. This is 
in agreement with the experiments because the gap in the diagonal direction is just equal 
to Ab2, introduced at T = T^. The minimum of |A(6')| on the diagonals is in agreement 
with other ARPES experiments@'i as well. One should of course notice that in our model 
closeness of two different A's only occurs in very narrow parts of the parameter space and 
therefore to a certain degree is accidental. 



VI. CONCLUSIONS 

We have formulated the one-loop renormalization group for a two-dimensional system of 
interacting electrons on a square lattice, described by the Hubbard model. For band filling 
different from one-half, the renormalization flow for / superior to some crossover /co comes 
only from the contribution due to the electron-electron bubble-diagram, while the electron- 
hole contribution decays exponentially as exp(— //2), where Ico depends on band filling, but is 
independent on the strength of the interaction. We decompose the BCS pairing interaction V 
for electrons in the vicinity of the Fermi surface in Fourier components of the five irreducible 
representations of the point group, defined at the Fermi surface. Diagonalizing V{lco) in 
each representation, we get five sets of decoupled BCS flow equations. The minimal (i.e. the 
most negative) eigenvalue of V{lco) determines the critical temperature and the eigenvector 
gives the form of the gap function. Unlike the usual approach,^! the characteristic of the 
procedure presented here to obtain the symmetry of the gap function is that only the axial 
coordinates at the Fermi surface is relevant, while the radial dependence is "scaled out". 
Moreover, the renormalization group treatment of the whole Brillouin zone, and not only 
of the narrow belt about the Fermi surface has allowed us to show that the origin of the 
attractive part of the pairing interaction in the Hubbard model is in electron-hole fluctuations 
on rather high energy scales, up to the bandwidth. 

We have calculated V{lco) for the case where the flow due to the electron-hole channel 
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can be treated perturbatively, i.e. when the filhng is far enough from one-half. The diagonal- 
ization of V in terms of angular harmonics gave us the type of superconducting instability: 
for weak filling, the instability occurs in the B2 (dxy) singlet channel, while for filling close 
to one-half, the Bi (dx^-y-i) singlet instability strongly overwhelms all others, what are the 
results in agreement with previous work.0@0 Particularly interesting is the fact that the 
order parameter Ab-^(^) can be very well approximated by the function cos2z{6), where the 
function z{6) depends only on the anisotropic Fermi velocity and on the geometry of the 
Fermi surface. This gives for Abi(^) a function that has peaks in the directions of the van 
Hove singularities. The slope of the peaks increases as we approach half-filling. This can be 
a justification to consider the interaction only between electrons in the close vicinity of the 
van Hove points as relevant if we are very close to half-fillingii. We believe that the form 
of Abi^{0) does not depend considerably on the strength of interaction, and recent T-matrix 
calculationsH for realistic values of the interaction Uo give a gap function in accord with our 
assumptions. 

We find a superconducting instability at any electron concentration away from half- 
filling. The underlying physical mechanism, namely exchange of spin or charge density 
fluctuations, is the same as in previous approaches. liHii We do however feel that our present 
results are on a more solid footing than the previous work because the present one-loop 
renormalization group scheme does not make any a priori assumptions about important or 
unimportant diagrams and provides a more systematic way of handling the dynamics of the 
fluctuations being exchanged. The only restrictions come from (i) the limitation to one- 
loop order, necessitating weak coupling, and (ii) the requirement that the e-h diagram are 
a perturbation with respect to the e-e diagrams, implying that we can not be too close to 
half-filling. The region of validity of the approach is shown in Fig. 5. Further, self-energy 
diagrams have been neglected, however, these are expected to produce important effects 
only at two-loop order, and therefore are expected to be negligible in weak coupling. 

In our weak coupling model the superconducting critical temperature is negligible (but it 
exists, for any filling!) if we are not in the immediate vicinity of half-filling, which means that 
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the Hubbard model with small Uo and small (perturbative) antiferromagnetic fluctuations 
does not suffice to describe the high critical temperature (~ 0.02t) of the cuprates. There 
exist two possible ways (related to two restrictions of our calculations) to increase the critical 
temperature. The first is to simply increase Uq and to remain far from the half-filling, 
keeping the e-h channel non-singular. To treat this case, an approach perturbative in Uq is 
possibly only of limited use. Ideally, renormalization should be done exactly, and not using 
a simple one- loop (or n-loop) scheme (which is actually just an "intelligent" version of the 
perturbative summation). We can speculate and suppose that even in the case of strong 
coupling there exists the crossover Ico, above which the fiow is of the BCS type. It is to expect 
that Ico is not very different from l^o that we have calculated. This means that the cutoff 
Tco for the effective BCS theory is Aoexp(— /co) (see Fig. 3. and eq.( |3.2| )), where Aq is the 
initial cutoff of the theory, equal to 8t. The second possibility to increase Tc is to approach 
half-filling very closely, making Fermi surface nesting important but remaining in a weak- 
coupling regime. For that simplified one-loop renormalization group calculationg22l 
has shown that superconductivity wins over antiferromagnetism only if the e-h contribution 
to the fiow decays before the divergence in the antiferromagnetic channel takes place. Thus, 
we can say that here too, the effective theory is of BCS type. The difference with the first 
scenario is that the effective cutoff is very small (Fig. 3), and that the coupling constant 
is very strong, due to the strong ffow in both e-e and e-h channels at all scales / < Ico- A 
very important feature of the nested case with a small Uq is that it can be treated in terms 
of the one-loop renormalization group, renormalizing U as a function only of three angular 
variables. This is allowed because all important physics (i.e. the majority of the e-e and e-h 
ffow) is contained in the vicinity of the Fermi surface, making the effective phase space to be 
a rather narrow square ring ±Aj; {8t ^ Aj ^ Aco) where the marginally relevant interaction 
is a function only of the angular position of the particles on the ring. 

Finally, we have discussed the possibility of a superconducting state with a mixed sym- 
metry. In the presence of Bi order, the ffow in the B2 channel (which is the second most 
attractive one for 0.206 < fx < 0.276) will be only weakly affected by a nonzero order param- 
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eter of Bi symmetry. This gives rise to two superconducting instabilities, with the critical 
temperatures for the Bi and for the B2 channel, and < T^.. We have given the form 
of the pairing function for the effective BCS theory. At T < T'^ the relative phase of two 
order parameters is = ±7r/2. The resulting form of the energy gap |A(^)| has no zeros and 
minima are in diagonal directions, providing a possible qualitative explanation of ARPES 
experiments by ShenS and the decrease of the anisotropy with decreasing TB 
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FIGURES 



FIG. 1. The one-particle irreducible diagrams for the vertices r2 and generating the 
renormalization of the self-energy and the interaction, respectively. 

FIG. 2. The e-h bubble (solid line) and its derivative over / (dashed line) for ii/At = —0.25 
(a) and /i/At — —0.02 (b), for momentum transfer qeh — 2ki^ || (tt, tt). 

FIG. 3. The scale Ix as a function of filling (n). For / > the e-h flow decays exponen- 
tially. The insert shows the relation between and (n). 

FIG. 4. The relation between the angular variable z and the observable polar angle 9 for 
-H/U = 2 X 10-"; n = 1(a), 2, 9(i). 

FIG. 5. The curve UoPehil oo) = 1. Below the curve, the e-h contribution to the 
renormalization can be treated perturbatively. 

FIG. 6. (a) The shape of the function Peh{l oo) in {z,z') space at —fi/At — 0.2. 
The nesting at half filling occurs for z = z' = ir/i. The split singular lines show the best 
incommensurate nesting vector, (b) The minimal eigenvalue of the pairing V{z, z') in every 
of 5 irreducible representations of D4 point group. 

FIG. 7. (a) The shape of three possible gap functions: Ab^(6') for /i/it = —0.001 (dot- 
dashed); Ab^{9) for fx/U = -0.5 (dashed) and Ae{9) for fx/U = -0.31 (solid line), (b) A 
very good approximation for Abi{9) is just cos22;(^), shown here for the same choice of /i 
as in fig. 4. The dashed line shows cos 9 to comparison. 

FIG. 8. The scale h — — InTc/St as a function of the logarithm of the chemical potential, 
for a few values of couphng Uq. For very small and for Uo/At > 0.5 the curves are out of 
the range of validity (see fig. 5). 
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